Geometric Laws of Vortex Quantum Tunneling 
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In the semiclassical domain the exponent of vortex quantum tunneling is dominated by a volume 
which is associated with the path the vortex line traces out during its escape from the metastable 
well. We explicitly show the influence of geometrical quantities on this volume by describing point 
vortex motion in the presence of an ellipse. It is argued that for the semiclassical description to hold 
the introduction of an additional geometric constraint, the distance of closest approach, is required. 
This constraint implies that the semiclassical description of vortex nucleation by tunneling at a 
boundary is in general not possible. Geometry dependence of the tunneling volume provides a 
means to verify experimental observation of vortex quantum tunneling in the superfluid Helium II. 



PACS numbers; 66.35.+a, 67.40.Vs 

The question how a singular Une with circulating flow 
around it enters a fluid in which friction is absent has 
pervaded the theory of superfluids for decades. It is cur- 
rently believed that at the absolute zero of temperature, 
a quantized vortex can come into the superfluid by a 
tunneling process. In this process, Galilean invariance 
has to be broken because there exists no absolute refer- 
ence frame in the homogenous bulk if the temperature 
equals zero. The formalism able to describe the quan- 
tum mechanical process of tunneling depends sensitively 
on the relevant length scales involved. If the scale after 
tunneling, when the vortex is able to move freely in the 
superfluid, is much larger than the core size of the vortex, 
we are able to resort to a hydrodynamic, i.e. semiclas- 
sical description jl]. If this scale is near the core size 
^, a description of vortex tunneling requires knowledge 
of the dynamic many-body structure of the vortex core, 
a formidable goal yet to be achieved. In this paper, we 
address the question of geometric laws involved in the 
exponent of vortex quantum tunneling in the semiclassi- 
cal approximation and investigate the limits imposed on 
this exponent by the shape of the flow obstacle breaking 
Galilean invariance. We show that the usual description 
of tunneling with a vortex moving in a pinning potential 
generated by some flow obstacle is limited by a geometric 
constraint related to the vortex core size. 

The tunneling rate of the vortex is proportional to 
exp[— 5'e(£')/?i], where Se{E) ^ ?i is the Euclidean ac- 
tion as a function of vortex energy E |0] : 



Se{E) = jdajdK ■ 



(1) 



Here, K = —iX{a) is the imaginary collective vortex co- 
ordinate and P(cr) the local vortex momentum per unit 
cr-length (with a the arc length parameter labeling points 
on the vortex string). The integral has to be performed 
along the classical vortex path with constant energy. We 
assume dissipation to be negligible in the very pure sys- 



tem under consideration. The only driving mechanism 
behind the tunneling transition we consider is the super- 
fluid Magnus force at zero temperature. 

We first show in its full generality that the action (|l|) 
in an incompressible fluid is a geometric quantity. The 
corrections induced by the mass of the vortex have been 
shown to be negligibly small in the semiclassical limit 
They will consequently be neglected here. 

In analogy to electrodynamics Q, the canonical mo- 
mentum P in an incompressible superfluid corresponds 
to the vector potential and is a gauge dependent vector 
with components 



Pr 



hpo b,jX'r 



(2) 



where hpQ is the bulk particle number density multi- 
plied by Planck's action quantum h, X'^ = dX^ /da 
the line tangent. We took the quantum of circulation 
to be K = {Ny/Ns)h/m, with Ny the winding num- 
ber of the vortex and Ns the number of real particles 
in the superfluid boson (a Cooper pair has Ng = 2). 
The antisymmetric gauge tensor b^i, of string dynam- 
ics 0, a 3-t-ld generalization of the stream function of 
classical hydrodynamics, is deflned by the dual transfor- 
mation v^ejafif = b^y,a + 6q^,i/ + Ka.fL from the con- 
servation of the four-current = pv^ (with p = po). 
Let ej, eg, Co- = X' be a righthanded co-ordinate ba- 
sis on the string. In our non-relativistic case, we have 
- §§ PadX'' hda = {Ny/Ns)hpa JJJy^dX^ AdX^ Ada , 
where g is the determinant of an arbitrary metric in the 
basis Ba, Ba [a = 1,2). This yields for the components 
of P in the directions , eg the relation 



d2Pi -diP2 ^ j^hpQy/g 



(3) 



Equation (|^) represents the conjugateness of positions 
and momenta in the directions ej , eg in its general form. 
As a well-known example, if we gauge Pg — 0, then, in 
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Cartesian co-ordinates, Px = {Ny/Ns)hpoY, the gauge 
to be used below in the application to the 2d case. We 
are now able to rewrite the action dll) as 



function of vortex position xi^Vi by i? = An E/mpQK^ , 
it takes the form 



S,{E) _ Ny 



(4) 



The new co-ordinate differentials are dZ^ = cosadK^ + 
sinadK^, dZ^ — — ainadX^ +coaadX^. The angle 
a{a) parameterizes the local (co-ordinate) gauge freedom 
for the momentum of rotations about Co-, contained in 
(^. In what follows, Ny/Ng = 1 for simplicity. The 
gauge invariant expression (^) tells us that the dimen- 
sionless action is the number density times the volume 
which is involved in the tunneling process. Note that the 
action is scaled with h, not h — h/2'K. This Volume Law 
is the dominant contribution to the tunneling exponent 
in a superfluid in which the non-relativistic limit is the 
physically relevant one. The number of particles in the 
effective tunneling volume plays the role of the quantum 
number in a Bohr-Sommerfeld type quantization of (|^). 
We stress that this result holds however complicated the 
motion of the vortex string actually is as long as we are 
in the semiclassical domain. 

What is of interest to us here are the limits which 
govern the actual value of the tunneling volume for a 
given external geometry. To tackle the general features 
of the problem, we restrict ourselves to the case of point 
vortices in 2d, where we can discard the complicated a- 
dependence in . We are then in the position to use the 
tools of conformal transformation to calculate the energy 
of the vortex as a function of its coordinates near some 
irregular boundary of the superfluid |^ . The simplest re- 
alistic geometry of interest to us for this purpose is that 
of a half-ellipse whose small semi-axis a is in the direction 
of the flow parallel to the boundary with a much bigger 
semi- axis & 3> a perpendicular to it. This provides a re- 
alistic pinning potential directly related to geometrical 
properties of the pinning site in contrast to a harmonic 
shape of the potential with fitting parameters used, e.g., 
in 

We can map the readily solvable problem of a vor- 
tex in the presence of a half-circle by the inverse of the 
Joukowski transformation [|j to that of the half-ellipse- 
vortex geometry in question. The vortex energy con- 
tains three terms. The first term is half the self en- 
ergy of a vortex pair and the second, which gives the 
pinning potential, stems from the two image vortices in- 
side the circle respectively ellipse required to satisfy the 
boundary conditions. The last one comes from a super- 
imposed external flow velocity u at infinity in the nega- 
tive X direction parallel to the boundary. We use appro- 
priate elliptic co-ordinates defined by x = Zsinhxcos?7, 
y = / cosh X sin ry, in which choosing x = Xo gives us 
a = I sinh xo , h = I cosh xo : with I = \J\p- — d} ~ b ^ S. 
the overall length scale. If we normalize the energy as a 
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a + b exp(xi - Xo)|sin??i| sinh(xi - Xo) 
^ (sinh^(Xi - Xo) + sin^ ViV^'^ 

smh(xi - Xo)| sm77i| . (5) 



We cannot solve for the path of constant E analytically if 
we allow for any values of ^, a, b and u in the semiclassi- 
cal domain. What we want to show here are general fea- 
tures with respect to the geometrical quantities involved 
in the tunneling exponent. Hence for the velocity u we 
consider the limit u <C n/iirl v^, where vl = 
is the characteristic velocity associated with the many- 
body quantum structure of the fluid ('Landau' velocity). 
The external current is limited by 2u/vl <C a/b, such 
that the velocity without vortex at the ellipse top re- 
mains well below the critical vl. 




FIG; 1. The two vortex paths of the same constant en- 
ergy Eq. The vortex on the path passing the eUipse top on a 
distance slightly larger than s > ^ has to tunnel to the sec- 
ond path, which is at a distance ~ Yo from the boundary. 
The small inset shows the shape of the barrier from (^ for 
rj — T]! — 7r/2. 



We now come to the crucial point of our argumentation: 
We require that the distance of the vortex to the el- 
lipse top, the point of closest approach of the vortex, 
is slightly more than f , because any smaller distance to 
the boundary would invalidate a description of the vortex 
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in semiclassical terms. This requirement can be imposed 
by Sxivi = ■'''/2) = s/a, with ^ < s <C a, and where we 
set Sx = Xi — Xo- The energy E = Eq oi the paths we 
consider (see Fig. |^) is then fixed to be 



Eq = E{a, b, s, ^, u) — hi 



h s 



47ru(a + h) s 

• (6 
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For low velocity, we can write expi?o — {b/a){s/£,), the 
value to be used in the following. 

There are two vortex paths having the same energy 
(Fig. |l|) . One of these is that of the vortex at the bound- 
ary and closely following the ellipse, for which Sx ^ I- 
We obtain from (||), if sin 771 sinh(5x, b/a 1: 



F^^^ 1 In 
Anu ox 



-Sx 

s 



(7) 



The other vortex path of constant energy Yn{xiiVi) is 
located far away from the ellipse and the boundary. The 
vortex has to tunnel through the barrier represented by 
(^ from the ellipse path (0) to this trajectory. In low- 
est order of the velocity u, this path follows a constant 
distance Yq to the boundary: 



47rw \ 27ru s 



(8) 



The two paths (|) and (0) do never meet provided that 
2u/vl <C a/6 holds. To describe the tunneling trajectory 
of constant energy, we have to use complex co-ordinates 
in the , space contained in (|^) such that they meet 
at some point in this space. We can then obtain a closed 
escape path. To the end of finding the path, we return to 
the original Cartesian co-ordinates and write the vortex 
trajectory in the complex {K = —iX, Y) space as 



fcosh^ixo + Sx) + tarLh '^{xo + Sx)K^ 



(9) 



where the deviation from the original ellipse Sx is in gen- 
eral via (|^) itself a function of K,Y. In lowest order 
of Sx/xoi however, (|^) yields approximately a constant 
Sx — s/a equal to its value at the ellipse top. The two 
trajectories (||) and (0) then cross at 



^ + 1)y,^PYo. 

b a/ 



(10) 



Evaluating the integral (||) leads to the tunneling action 
(see also Fig. ||): 



= 2po (^N - Y^)dK ~ poI3YS 

K a/b 
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27ru s 



(11) 



We remind the reader that the result in this simple ana- 
lytical form is valid to lowest order in the small quantities 



2u/vl ^ a/h <C 1, s/a <C a/h. But the feature that s, a, b 
enter the action does also hold for more general values in 
the parameter space. We conclude that the semiclassi- 
cal action depends on three geometrical quantities: The 
characteristic lengths parallel and perpendicular to the 
flow a, b and the closest approach distance s to the ellipse. 
The quantity s enters because we imposed the condition 
that the vortex must not get nearer to the surface than 
^ < s. This resulted in an effective rescaling ^ £, exp Eq 
of the ultraviolet cut-off in the vortex energy logarithm 
and a restriction on the effective sharpness /3 > a/b-\-^/a 
of the ellipse. 




FIG. 2. The area ~ I3Y^ = "Kj^^^ enclosed by the 
vortex path in complex space gives the action (^ij) in two 
space dimensions. This area has a lower semiclassical limit 
{a/b + aa)Yo^- 



As a generalization of (|l l|) , the semiclassical action for 
a massless vortex in d dimensions will take the form 



(12) 



where /? ^ 1 characterizes the effective dimension of the 
vortex escape path, i.e. the relative degree to which it 
is confined to n dimensions by the presence of an asper- 
ity which is effectively n-dimensional. The quantity /? is 
bounded from below by the requirement of semiclassical- 
ity for the vortex path at the boundary, as expounded 
above for the two-dimensional case. The tunneling vol- 
ume vj^f'^ is that for a vortex escape path of 0(d — 1) 
symmetry pO[ |, which is the highest possible symmetry 
if one preferred direction, namely that of the external 
current, is given. In our d = 2 case of (pdl), we have 
n=l, /?~a/5-f s/a and = Yq. In three dimen- 
sions {d = 3) we have V^^^ = (47r/3)i?o for ^ half-ring 
with radius Rq. 

The conclusions we have drawn should be able to shed 
light on the question if in experiments on the flow of 



3 



superfluid 'Hie through sub-micron orifices we are deal- 
ing with the quantum tunnehng of vortices below about 
T 



150 mK 1 12 
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Geometry dependence provides a 
tool to verify that the critical velocity plateau observed 
there comes indeed from a quantum process and not an 
alternative process of classical instability, e.g. of the 
Kadomtsev-Petviashvili type [Q, is taking place. In "fie, 
the velocity wl is quite large due to the small coherence 
length in the order of the interparticle spacing (it equals 
the Landau velocity « 59 m/sec at p — I atm if ~ 2.7 K 
is taken) . The observed critical velocities of flow through 
the orifice on the temperature independent plateau are of 
the order of 10 m/sec [13 1^. This leads to materializa- 



tion radii Rq of the order of nanometers {Rq corresponds 
to the 2d quantity Yq in (H))- A hydrodynamic treatment 
appears thus justified and there should be a strong de- 
pendence of critical velocities on the surface structure of 
the orifices. The simple model of vortex half rings situ- 
ated at plane walls with their momentum axis antiparallel 
to the flow (l^Jl^ does not account for the required ex- 
plicit incorporation of Galilean invariance violation. The 
enlargement of the core size at boundaries used in this 
model to explain observed critical velocities is naturally 
contained here. It is represented by the condition that 
a vortex moving in a pinning potential caused by a non- 
trivial geometry only exists as a semiclassical vortex able 
to tunnel through the barrier if its distance to the wall 
is everywhere bigger than ^. This results in a nonzero 
tunneling energy i?o > of the vortex interpretable as 
an enlargement of ^ ^ ^ cxp Eq . 

The necessity of nonzero tunneling energy is responsi- 
ble for the fact that we do not deal with the nucleation of 
vortices here, i.e. their creation from 'nothing'. The vor- 
tex already has to exist at the boundary with a small but 
nonzero energy such that it can tunnel to get a free vor- 
tex. Hence we are dealing with a depinning rather than 
a nucleation event. In the semiclassical limit the descrip- 
tion of the nucleation of zero energy vortices is in general 
only possible if the curvature radius of the flow obstacle is 
constant, that is in the case of d— 1-dimensional spheres. 
If the curvature radius of the obstacle varies, the distance 
of the zero energy vortex to the boundary varies and can 
get smaller than ^ at the point of closest approach, thus 
violating the requirements of the semiclassical approach. 

To conclude, we emphasize again the important fact 
that violation of Galilean invariance in the direction of 
the flow is a necessary condition for vortex quantum tun- 
neling to happen at the absolute zero of temperature. 
The tunneling exponent will depend on geometrical prop- 
erties of the vortex path associated with curvature scales 
of the asperities breaking the invariance. That the semi- 
classical description be valid introduces the additional 
length scale s, the closest approach distance to the flow 
obstacle. The semiclassical theory of vortex quantum 
tunneling in a pure superfluid is essentially a geometric 
theory. 
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